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Abstract 

The pattern of spontaneous chiral symmetry breaking (CSB) in confining back¬ 
ground fields is analyzed. It is explicitly demonstrated how to get the inverse square 
root large proper time asymptotic of the operator product expansion which is needed 
for CSB. 


1 Introduction 


The phenomenon of spontaneous chiral symmetry breaking (CSB) is one of the two 
most important nonperturbative properties of quantum chromodynamics (QCD). An¬ 
other property, confinement of color, is believed to be deeply connected with CSB. It 
is interesting to see how conhnement and CSB sometimes play competitive roles when 
one speaks about low-energy hadron physics. Namely, one can hnd statements in the 
literature that confinement is ’’not seen” in low-energy spectroscopy, and it is possible 
to model the QCD vacuum by a set of classical field configurations (as one successfully 
does, for example, in the instanton model P, see reviews 1210 ] and references therein), 
despite the fact that such ensembles typically have no confinement in the sense of area 
law of the Wilson loop. There is also an alternative line of arguments, going back to the 
good old constituent models, which states the ’’supremacy” of conhnement (understood 
in terms of the conhning string formation and corresponding linear potential between 
quarks), while CSB is to be derived from conhnement (see, e.g. j3]) and there is no need 
in instantons or alike in this approach. In fact, these two points of view are to some extent 
complementary to each other, as the example of pion clearly shows: pion is simultaneously 
Nambu-Goldstone boson which should be massless in the chiral limit and, on the other 
hand, it is a bound state or quarks and antiquarks, as any other meson in QCD is. 

The two physical pictures outlined above correspond to two diherent approaches 
one usually uses analyzing the phenomenon of CSB in QCD. The hrst approach concen¬ 
trates on studies of the Dirac operator zero modes in this or that gauge field background. 
The chiral condensate is related in this case to the (quasi)-zero modes density by the 
well known Banks-Casher relation [Sj. The model dependence enters when one starts 
to answer an important question about physical relevance of the chosen background, in 
particular, about its conhning properties. We refer the interested reader to the review [2] 
and references therein for the details of this approach. 

Alternatively, one can choose this or that interaction kernel (for example, manifestly 
providing conhnement) and construct gap equation for the chiral condensate with this 
kernel (see, e.g. recent paper jB]). The same procedure is commonly used in hnite density 
QCD [7j. This exhibits in a very clear (but mo del-dependent) way the general relation 
between conhnement and CSB jH], taking conhnement as the cause of CSB. 

In Shifman-Vainshtein-Zakharov sum rules approach j2] the chiral condensate {qq) 
enters as an input parameter not directly related to, for example, gluon condensate {agF"^). 
In other words, there is no simple relation like {qq) = const ■ {asF'^)^F one might naively 
think of.^ In fact, studying correlators of hadronic currents one can get more sophisticated 
relations (see, e.g. HDl ) between chiral and gluon condensates (and other nonperturbative 
quantities like /j^). The problem however is that SVZ approach misses the relation between 

^It is worth reminding that nonperturbative gluon condensate is not a local order parameter, while 
chiral condensate is. 
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chiral symmetry breaking and confinement. For example, nonzero value of (agF^) does 
not at all indicate that the vacuum confines.^ As is well known, one possible way to take 
the effects of confinement into account properly is to consider the dynamics of nonlocal 
objects like the Wilson loops. 

In the present paper we address the problem of CSB in the first-quantized language, 
i.e. in terms of quark trajectories and not fields. In this sense, we are closer to the 
latter approach discussed above and not to the former one. The criterium for CSB in this 
framework is given by well known Banks-Casher asymptotic law (eq. Q of this paper). 
We are going to address the following question: how this asymptotic law follows from 
the properties of the Wilson loop expansion over local condensates and their derivatives. 
Speaking differently, we are looking for the simplest subseries of the Wilson loop operator 
product expansion, whose summation provides CSB at the level of one-loop effective 
action. It is shown that CSB has to do with the large proper time asymptotic limit of the 
specific confining nonlocal gauge-invariant correlator of gluon fields. 


2 One-loop effective action 


We start from the expression for the Euclidean QCD partition function: 

Z = J VA^VqVq^ exp J + j dW - im) q^ = 

= {det{i'y^D^-im)) (1) 

and we confine our attention to the case when the mass matrix is proportional to the unit 
matrix in flavor space with the eigenvalue m. As usual, the average of any gauge-invariant 
operator 7^(A) over gauge fields is understood with the standard Yang-Mills action 


(7^(A)) = 


7^(A) exp 



d'^xF, 





( 2 ) 


where normalization factor, proper gauge-fixing and ghost terms are included in the inte¬ 
gration measure. The chiral condensate is given by the standard form 


(gg)W) = ^{q^q)F) 


1 dlogZ 
V dm 


dTeff 

dm 


( 3 ) 


and the superscripts M, E stays for Minkowski and Euclidean values. It is assumed that 
the right hand side of © does not vanish in the limit m —> 0. This nonzero condensate 
corresponds to spontaneous CSB. 


We are going to exploit the so-called Feynman-Schwinger representation technique, 
whose essence goes back to the seminal papers nacni. There are basically two lines of 

^As one can see on the lattice at large temperatures m 
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use of this approach in the modern research. The hrst one exploits advantages the path 
integration provides for gauge-invariant formulation of relativistic bound state problems 
(see review ^3] and references therein). The second line |151 ITT)] concentrates on the 
loop calculations in perturbative held theory (for review see |I7j), which is sometimes 
simpler in world-line approach than in conventional Feynman diagrammatic framework. 
We collect some relevant formulas in Appendix A for convenience, while the interested 
reader is referred to the original papers, textbooks UBEHEII and cited reviews for all 
technical details. 

Conhning ourselves by quenched approximation (exact in large Nc limit), we have 
(det iF) = exp(logdet iF) -|- 0{N~^) and hence the standard expression for Euclidean 
one-loop effective action: 


1 /•oo ^7^ 

Te// = — logZ = -2 —exp{-m‘^T) ■ {Z[A,T]) (4) 

where 

ZlA, T] = J VXf, J exp (-S'o) tr P exp J dr (A^i^ - ?/’,.) j (5) 

with the free world-line action given by S'o = / dr • We have included the 

factor \jV in the integration measure 'Dx^. The factor 2 = 4 x (1/2) in came from 
the trace over anticommuting coordinates (i.e. integration over the fermionic helds with 
the free action is normalized to unity in our conventions). All dynamical information 
is contained in the double average (over gauge helds and over quark trajectories) of the 
spinor Wilson loop 

u;(T)^(Z[A,T]) (6) 

There are several limiting cases where one can successfully study the behavior of w{T) or 
related functions. Of prime importance is the nonrelativistic limit of large mass m, where 
the target space for the contours x^ becomes ehectively three-dimensional and introducing 
the einbein helds for dynamical masses one can systematically explore constituent picture 
of hadrons (see, e.g. nain] and references therein). However for the trace (and hence 
for the closed contours) we are discussing at the moment nonrelativistically suppressed 
backward-in-time trajectories are as important as forward-in-time ones and to address this 
problem in einbein helds formalism one probably has to use some alternative methods (see, 
e.g. m)- 

The second important case corresponds to the small T asymptotic. In a way, this 
is the standard operator product expansion [22] • In context of the theory of gravity it 
corresponds to well known Schwinger - DeWitt expansion m One is to expand the 
Wilson loop in powers of helds in this limit. For constant background helds this is the 
way one obtains the ehective Lagrangians of Heisenberg-Euler type. The typical term of 
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this expansion looks like 


l^jjkx pm\£)k2 pm2 pkp prup^ _ pi 


( 7 ) 


The leading term represents the so called heavy quark condensate jH] 

(see also where the heavy quark condensate is discussed in the path integral formal¬ 
ism, including nonzero temperature case). It is clear that this expansion is unapplicable 
in the limit of vanishing mass m. 


The phenomenon of CSB is related to the large proper time asymptotic of w{T). As 
it was noticed by Banks and Casher in their seminal paper [^, for spontaneous CSB one 
should have 




at T —cx) 


( 9 ) 


Indeed, it is easy to see from and 0 that 


(qq) 


dTeff 

dm 


rsj 



dT exp(—m^T) 


c 

vf 


—d-y/vr • c 


( 10 ) 


Taking into account that in free case (i.e. without gauge fields) w(T) is given by 


wo{T) = (47rT)-i 


( 11 ) 


one can say that the quark is dynamically forced to move effectively in 1 dimension instead 
of 3 -|- 1 and this is the cause for CSB in this framework. 


There are a few well known cases where the dimensional reduction of this kind indeed 
takes place. However, it is worth noticing that the condition m is rather restrictive to be 
incorporated in a simple way into the standard background held formalism. Indeed, one 
can show on general grounds (see |2fil 1771 |2H] and [22] for review) that the typical large T 
asymptotic of the heat kernel trace Tr/C(T) = / dx)C{T, x, x), entering the effective action 
as 

CX) 

= ^ / '=(r) (12) 

0 

in generic fixed background held is given by 

where nonlocal factors Wn can be expressed as integrals of the corresponding zero modes. 
Notice that the expansion goes in integer powers of 1/T. We will address this contradiction 
between m and 0 in the Section 4. 
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Simple illustrative example is Heisenberg - Euler effective Lagrangian ra for constant 
magnetic field 

CO 

^ = -^/ ^exp(-m^T) ■ cth eHT - ^ (14) 

0 


The lowest energy level of the massless fermion from Dirac sea in constant magnetic field 
is zero, hence the absence of exponential damping in (I14|l . cth eifT —> 1 with T —oo. 
The factor 1/T = T ■ in front of the cth eHT corresponds to the fact that in 

four-dimensional space-time there are two directions the fermion can move along as a 
free particle, while the dynamics in two other directions is confined by the field.^ This 
is nothing but the first term in the rhs of m- On the other hand, we see something 
new here with respect to m- The leading term at large T is the last term in the rhs 
of ()14|1 . which is 0{T^) and not 0{T~^). This term represents one-loop short distance 
charge renormalization and, at the same time, the leading strong-held (i.e. small mass) 
logarithmic asymptotic of (HI: 


L 


24 ^ 


(15) 


This phenomenon of strong held - short distance duality (small T - large T duality in 
our context) is quite general (see, e.g. recent discussion in [S2]) and provides interesting 
possibilities for OPE subseries summation (see discussion below). Thus we see that quan¬ 
tum dynamics (the necessity to express the answer in terms of renormalized quantities in 
this case) can make the result m inapplicable. 


3 Effective action at Gaussian level 


The spinor Wilson loop factor Z[A, T] given by (0) can be expanded in powers of helds and 
derivatives. It is convenient to use Fock-Schwinger gauge condition, which is a particular 
case of the so called generalized contour (or coordinate) gauge [221 • The latter is dehned 
in terms of the oriented non-selhntersecting contour z^{s) as 

( 16 ) 

The simplest contour gauge one usually uses is the Fock-Schwinger gauge condition with 
Zfj_{s) = -|-s(x —In this gauge d'^z^jds^ = 0 and thanks to that one can express 

the vector-potential in the following form (compare with (I58|l from Appendix B) 

1 

= / s^s2/pexp(s?/^D^)“^F^^^(x(°)) (17) 

_o_ 

^For the field H = e^iJ, the particle moves freely in z-directions (along the field) and in t-direction 
(no force acts on the particle in rest). 
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where y = x — and Latin indices a,b = 1 ,— 1 stay for adjoint color. 

It is worth stressing that for contour gauges gauge invariance corresponds to the con¬ 
tour independence, not just independence, as is sometimes posed. This situation is 
analogous to the covariant i?^-gauges, where ^-independence is necessary but not sufficient 
condition of gauge invariance; in other words one can easily construct gauge-noninvariant 
(and hence physically unobservable) but ^-independent quantity. The contour indepen¬ 
dence is restored only in the full Wilson loop, but not at any given order of the expansion 
over fields and/or derivatives (neither at any given order of the covariant perturbation 
theory 123 EH). Practically it means that general expansion of some nonlocal object like 
the Wilson loop over local condensates has no universal coefficients, independent on the 
choice of contours used to fix the gauge (see the Appendix B). 

With these reservations in mind, we can proceed and expand the Wilson loop over 
gluon fields. In quantum case the dynamics is determined by the average over quark 
trajectories 2 /^(t) (and over spinor ’’coordinates” with Gaussian weight exp(—A q) 

and over the vacuum gluon fields with the standard Yang-Mills action. Having performed 
the latter average, the first two nontrivial terms in the expansion of Z[A, T] take the form 

T n 1 1 

{Z[A,T]) ^ wo{T) + J Vy^ J exp(-S'o) J dn j dT 2 j dsi J ds 2 

0 0 0 0 

{siy^^{Ti)yp^{Ti) - - Si)) ■ {S2yp^{T2)y p^{t2) - - S2))- 

■ (Tr exp(s 2 l/(T- 2 )T>)iV 2 P 2 (a^^°^)exp(si|/(ri)T))F^jp,(x(°^)) (18) 

The computational technique for such integrals is well developed [THl [THl IT71 122] . The 
basic ingredients are the one-dimensional Green’s functions on a circle Gb{ti,T 2 ) and 
Gf{ti,T 2 ) defined by 

{yp{Ti)y^{T2))y = -Sp^GB{n,T2) = ^|ri - T2\ - ^ (19) 

and 

{'ipp{ri)'4)y{T2))i, = ^GF{ri,T2) = ^sign(ri - T 2 ) (20) 

The following identity is of special use 

(exp {v,(T,)km) exp {y,{T^)k‘P) = exp (-Gb(t„ ) (21) 

The result is given by the expressions 

/*00 

hh = i, ^exp(-m2r)TO„(r) . r-K(T) (22) 

where Wq{T) is defined by (fTT]l . while K{T) reads 

K{T) = (Tr (23) 
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with the formfactor 


i 

JF(^) = J duu{l — u) exp{u{l — u)^) 


(24) 


and ^ = TD„D„. The arrows indicate that the derivatives act on the right. The contrac¬ 
tion of indices in (ESI) is worth noticing (see ESI in this respect). The most nontrivial thing 
is the exact expression for the formfactor (I24|l . which was obtained for the hrst time in 
PHj for classical backgronnds (and effectively reprodnced in jSSj nsing Feynman-Schwinger 
techniqne). For large m eqs. lead to the expression (jH)). 


4 Large-T asymptotic limit and CSB 


We are to stndy large-T asymptotic of (1^ . One has ~ 1/^^ as —*• —oo. It is 
easy to show that for the scalar particle one wonld have ~ 1/^ as.^ —*• —cxo. This 
important difference corresponds to the fact that no exponential damping at large T other 
than exp(—m^T) is possible for fermions in the chiral limit. Naively one has 

lim T^KiT) = 2(Tr F^,(xW)D-"F^,(a;(°))) (25) 

T—^oo 

This expression is formal, however, dne to infrared divergencies. We will show that general 
asymptotic expansion at large T may contain logarithmic terms of the form 

T27F(T) = colog(J)+0(T-MogT) forT^cx) (26) 

where A is typical correlation length and Cq - some dimensionless coefficient. Possible 
snbleading T-independent contribntion in the left hand side of (I26j) is inclnded into the 
dehnition of A. 

Let ns show how the behavior (EHI) of K{T) follows. The asymptotic pattern is 
controlled by the fnnction 

/(s) = (Tr F^.(a;®) exp (27) 

which one needs to know both at large and at small s. The small proper time asymptotic 
of /(s) is given by the standard OPE: 

f{s) = (Tr Fl,) + s(Tr F,,D^F,,) + ^(s^) (28) 

As usnal in SVZ snm rnles, it is assnmed by dehnition that all pertnrbative contribntions 
are snbtracted from each term of EHl), thns dehning the gennine nonpertnrbative fnnction 
EB- In the framework of Wilson OPE each term in EHl) depends on the dynamical 
scale (separating contribntions of pertnrbative coefficient fnnctions and nonpertnrbative 
matrix elements). The snbtle qnestion abont /i - dependence of the fnnction /(s) and 


the effective action is somewhat beyond our main line and will not be discussed. Taking 
more phenomenological attitude, the reader may think of f{s) as being computed on 
nonperturbative field configurations of one’s favorite QCD vacuum ensemble (instantons, 
dyons, P-vortices etc); this would correspond to some effective /i ~ 1 GeV of the order of 
the onset of nonperturbative dynamics. 


The Gaussian approximation to the function /(s) defined by (1771) can be written as 
(see details in Appendix B) 


/(s) ^ /2(s) = 


(dvrs) 


dHexp {-f/4s) (Tr F^,(a;(°))exp (29) 


The correlator in the right hand side (compare with (ITTIl l is frequently used in the Gaus¬ 
sian stochastic scenario of confinement ((SH], see review [25 and references therein). At 
distances I larger than some typical correlation length A this correlator decays exponen¬ 
tially^ with 1. Therefore for sX~^ 3> 1 the leading asymptotic is given by /(s) ~ To 
be more precise, it is convenient to parameterize two-point correlator in the standard way 


(Tr F^,(a;(°)) exp Fp.(a;®)) = + 




- 


{lpSp.-Upp)D,{f) 


(30) 


where both functions and Di{P) exponentially decay at large distances. Gorre- 

spondingly, one has 

/2(s) = /d(s) + /di(s) (31) 

with the following asymptotic limits at large s: 


/dW = (5 + 0(s-’) ; /D,(s) = i + 0(s-‘) 

The nonperturbative constants t], ( are given in this approximation by 


CXJ 00 

= ; ( = - j dl^ 1‘DiU^) 


(32) 


(33) 


The difference between asymptotic limits of /z)(s) and /dj (s) is of crucial importance. 
Indeed, combining (17^ . and (IH^ one gets 

1 

K{T) = KoiT) + KdAT) = jdu u{l - m)[/d(m(1 - u)T) + /b,(m(1 - u)T)] (34) 

0 

and for T —> cx) we have 

T^Kd{T) ~ logT , while T^KdAT) ~ ^ (35) 

"^Notice that this length A characterizes the function D{1^) and, generally speaking, A 7 ^ A. On the 
other hand, it is physically natural to assume that A and A are of the same order of magnitude. 
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The latter term being exponentiated cannot produce term and hence the function 

alone gives no spontaneous CSB. In other words, D{1^) = 0 implies unbroken chiral 
symmetry. On the other hand, it is well known [211 that just nonzero D{t^) is responsible 
for confinement, while Di{P) is not. 

It is reasonable to expect that the pattern we have discussed at the level of two-point 
correlator is general. From higher correlators one would have ~ logT terms, which add 
up to the coefficient r], and ~ log” T terms coming from the corresponding reducible part 
(i.e. the product of lower order correlators). The exponentiation of the series produces 
the desired power-like behavior of w{T). All terms ~ T~^ logT are subleading and do not 
change the leading power. 


Of course, we cannot compute from the first principles the coefficient in front of a 
general log” T term of this series to make the above arguments quantitative. But confining 
ourself by Gaussian approximation we can proceed further. The actual results for the 
chiral condensate and other quantities will crucially depend on the profile of the function 
/ 2 (s), which, in its turn, is determined by the functions and Di(/^). Namely, from 

(P7|) and PI|) we have 


T^Ko{T) = - 


dlVD{l‘ 


du 




u{l — u) 


exp 


4Tm(1 -m)^ 




= - / exp iFo = 


2T 


2T 


3 

2 


jilVDiP) log {§,) + A + o(T-,ogT) 


Thus the leading large-T asymptotic has the form 

T^Ko{T) ^2r^\og(^^ 
with the correlation length A defined by 

OO 

f dPPD{P)log(e~'P/4) 

logA" = i 


/ dPPD{P) 


(36) 

(37) 

(38) 


It is worth repeating that this leading logarithmic term is absent in the deconfinement 
phase where D{p) = 0. 


To summarize, the leading large-T asymptotics (EHl) of Te// in conhning background 
from the Gaussian term reads 

^ exp{-m‘^T)wo{T) (^1 + 2rilog + ...^ ( 39 ) 
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where both parameters rj and A are of essentially nonpertnrbative origin. 


It is again instructive to compare dSi with the constant field case. For constant field 
the (unrenormalized) result (jl4|l can be written as 


r = -2 


dT / 

— exp(-m^T)wo(T)exp ( ^ Kn{eHT) 


2 n 


\n=l 


where 


K, 


= 2 


(-1)-+! C(2n) 


r2n 




(40) 


(41) 


n TT^ 

The ’’correlation length” for the constant held is inhnite and in this sense T is never 
large, the series in n cannot be truncated and at the Gaussian level K{T) —>■ const ^ 0 
for T ^ oo. For hnite correlation length typical term in the corresponding expansion 
does not increase as a power of T and K{T) ~ logT/T^ —>■ 0 as T —cx). In terms of 
PI) the duality jHI] mentioned above corresponds to the fact that the leading large T 
asymptotic of renormalized effective action (i.e. the expression m) is controlled by the 
lowest term of the unrenormalized expression (HOD (i.e. the coefficient ki). 


The leading Gaussian large T behavior in conhning background given by P|) is much 
softer that constant held answer (PD. It is clear that the logarithmic term (PI) as it is 
cannot lead to GSB (since we need power-like rise to get (jH)))- But partial summation of 
such terms can do the job. The most crucial point is the structure of the series in the 
right hand side of dSi. No universal closed form expressions analogous to (P|) . (P|) for 
all terms of higher orders are known (see, however, m for explicit form-factors of the 
third order). On the other hand, one can argue on physical grounds that terms 



(together with other ones) should present. Such terms correspond to the factorized part 
contribution of the higher order averages {Hi FF...F). The important role played by these 
factorized averages is known under the name of vacuum dominance for a long time and 
it is successfully used in sum rules approach. From general analysis of [S] we expect that 
it is the confinement property that causes the GSB and it is known for a long time that 
just described reduction of the gluon ensemble (known as Gaussian approximation in the 
context) provides conhnement (see review |2I] and references therein). Therefore one can 
hope that we have summed ’’many enough” terms to keep GSB. On the other hand, the 
connected parts we have omitted physically correspond to the exchanges by multi-gluon 
glueballs and gluelumps, which are heavy objects and hence their contribution is to be 
suppressed for the low energy physics (see Ea in this respect). Another argument comes 
from the abelian dominance picture (see review j3H| and references therein). Since in the 
maximal abelian gauge the higher irreducible correlators are suppressed (because they 
correspond to diagonal - off-diagonal gluon couplings), one gets the same factorization 
pattern. It is important that the Gaussian factorization is to be assumed for the integra¬ 
tion over as well (i.e. we speak about some kind of ’’rainbow” approximation and no 
contraction of belonging to different clusters is done). 
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Summing of this ’’Gaussian” subseries of the full Wilson loop average would result 
in the effective action 


roo ^ 7 ^ 

^eff = -2 — exp{-m^T)woiT) 


exp{T^K{T)) + ... 


(43) 


where dots stay for the contributions of non-Gaussian terms. The spontaneous chiral sym¬ 
metry breaking condition has the following form in considered Gaussian approximation: 


lim T-^T^K(T) = -—^ (44) 

T^oo dT ^ ’ 2 ^ ’ 

or, in terms of (El 

OO 

J df fD{f) = 1 (45) 

0 

in four dimensional space-time. Then for the condensate one gets 


(w> 




( 46 ) 


where A is dehned by (El- The condensate vanishes in the deconhnement phase transition 
point. 


The above result deserves a few comments. Physically, the parameter rj is given 
by some integral moment of the function D{P) and, at the hrst look can be of arbitrary 
value, while we need strictly ?7 = 3/4 to get GSB. In a sense, this is an artefact of Gaussian 
approximation. It is worth repeating that higher terms in the r.h.s. of El bring both 
~ log T terms, shifting pure Gaussian value of t], given by El to its pure ’’geometrical” 
value El and ~ log"'T terms, adding up to the exponent. To some extent it resembles 
well-known effect of unphysical surface dependence of the Wilson loop average computed 
in Gaussian approximation: to get rid of it one has to sum the full cluster expansion.^ 
Moreover, suppose that the Gaussian asymptotic law for dH would be different from 
log T law El- The conclusion in this case would be that Gaussian approximation does 
not provide GSB and one needs the full series (or some subseries other than Gaussian) 
in dH to get dl- The actual conclusion is different: Gaussian reduction of the conhning 
vacuum can lead to GSB (because of (El) if it is done in a self-consistent (in the sense 
of dm), dH) way. 

On the other hand, all dimensionfull quantities in nonperturbative theory (like con¬ 
densates, string tension, correlation lengths etc) are proportional to the corresponding 
power of Aqcd with some dimensionless coefficient, unequivocally hxed by the theory. 
Neither the ratio of such coefficients is a freely adjustable parameter. Since we discuss in 
this paper not the full theory in the gluon sector but rather its Gaussian reduction (still 

®This analogy should not be understood literally, of course, there is no any special surface in the 
discussed problem, since we sum over all trajectories of the light quark. 
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keeping confinement), the relation ()45|1 can be understood as a kind of self-consistency 
condition for Gaussian approximation (if we wish it to provide CSB). The fact that such 
self-consistent reduction should introduce relations between condensates and correlation 
lengths is well known The ultimate reason for that are the Bianchi identities: correla¬ 
tor of derivatives with any operator, for example, Fafsiy, (inversely 

proportional to some typical correlation length) is given by the higher correlators of F’s 
(another term, containing Bianchi form, vanishes). But the higher correlators factorize 
to the product of Gaussian ones in the chosen approximation. 

It is clear from the discussion that we were mostly interested in qualitative pattern 
of the effect. Nevertheless it may be interesting to compare the results with the exist¬ 
ing phenomenology of Gaussian stochastic picture of QGD vacuum. Standard way of 
representing the lattice data on the nonperturbative function D{F) is 

D{1^) = D{0)P (l/X) exp (-//a) (47) 

where P{x) is some rational polynomial, normalized by condition P{0) = 1. Lattice ^Tj 
data correspond to A = 0.22 Fm for quenched SU{3) case, while the value of dimensionless 
product iA(0)A^ is numerically about (0.15 ± 0.05). Unfortunately, the poor knowledge 
of the pre-exponential factor P{x) precludes one to make quantitative predictions from 
the expressions (USD, dH). However, it is interesting to notice that the simplest ansatz 
P{x) = 1 does not describe data satisfactory;® it predicts iA(0)A^ = 1/12 and too small 
value of the chiral condensate {qq) = —(140 MeV)^ instead of the correct value {qq) = 
— (250 MeV)^. Of course, one can easily £t P{x) to get the desired numbers. But this 
seems to be misleading since, as it has been already mentioned, eventually these are the 
higher order terms which shift r] and {qq) to their correct values. 


5 Conclusion 


We have studied the chiral symmetry breaking at one-loop level in the background of 
conhning gluon helds. The latter is taken in the Gaussian picture, i.e. we have omitted 
higher than two-point irreducible condensates. This approximation is supported by the 
sum-rule phenomenology as well as by more sophisticated analytical and lattice analysis. 
It is worth stressing that the vacuum gluon helds ensemble reduced in such way still has the 
conhnement property. It is shown that this vacuum breaks chiral symmetry spontaneously 
provided the large proper-time asymptotic of the operator product expansion has the form 

P|). 


At the same time it is worth stressing that despite ’’Gaussian reduction” alone is 
enough to get GSB, it does not mean that higher order terms are not important. There 

^Notice that this ansatz approaches the point I = 0 linearly: D{p) — D{0) ^ -s/F instead of the 
behavior D{P) — D{0) ^ P dictated by it^ 
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are at least two respects where they are: hrst, their presence is a matter of principle 
for correctness of the general Wilson loop asymptotic (see, e.g. comments in j2Ij) like 
snrface independence etc; and second, they can be numerically signihcant (say, contribute 
10-15%) in the quantities like r], (. 

We have also discussed that perturbative asymptotic of the corresponding correlators 
cannot provide CSB; roughly speaking, at small momenta perturbative Green’s functions 
are ’’too soft”. No any specihc topological properties of the vacuum gluon helds show up 
in our analysis, the only relevant thing is the large proper time asymptotic of the OPE 
(resulting from conhnement). We did not use in any place explicit expressions for the 
gluon helds prohle, only the vacuum averages entered our analysis. Nevertheless it would 
be interesting to understand the relation of the discussed issues with the well established 
phenomenology of CSB in, e.g. instanton backgrounds (without conhnement per se), 
having in mind that it is ultimately the strong conhning color forces that determine the 
dynamics of CSB in accordance with general analysis of jH]. 
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Note added 

The related issues are discussed in the recent preprint hep-th/0511176 by V.Vyas we 
got after our paper had been completed. The main interest of the cited paper concentrates 
on the consequences the asymptotic law © has for the behavior of two-point correlators of 
hadronic currents, while we have tried to answer another and physically different question: 
how can the asymptotic Q itself be obtained using the OPE language. 
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Appendix A 

Throughout the paper we use Feynman-Schwinger world-line integration technique. 
The standard building blocks here are the following: 

1 ) 75 - invariance of the determinant (we dehne — iA^)\ 

logdet + "i) = ^ logdet + m] + m]) = 

= ^ log det [ 7 ^ 7 ^] + (48) 

2 ) proper-time representation of the logarithm; 

7 dT 

\ogA=-J —exp (-AT) (49) 

0 

The above expression is of course symbolical and properly regularized form of (gni) must 
be used in actual computations. The typical examples are given by exact equality 

00 

log A -I- 7 = — lirn J dTT^~^ (1 log T) exp {—AT) (50) 

0 

which can be obtained differentiating the identity A“^r(l + 0 = /o°° dTT^~^ exp(—AT) 

or by frequently used integral 

7 dT 

logA =-y —(exp(-AT)-exp(-T)) (51) 

0 


3) integration over commuting paths replacing the bosonic traces: 

f d‘^p 


Tr exp [—T{—D‘ 


= / To: exp 


(2vr) 


(p|exp (-(p +c/A)2t) Ip) 


J I trPexp iig J A^(x)i:^dr 


(52) 


4) integration over anti-commuting helds representing the path-ordered exponent of the 
gamma matrices: 


Tr P exp 


' 4 


77„| I = / DV’trP exp I - / <ir ( -Vv’/’e “ 


(63) 


where the small trace tr goes over color indices only. The corresponding operators i/j 
anti-commute and act on the Hilbert space of Dirac spinors ja) as 






( 54 ) 
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The normalization is provided by 


J T’V’exp j dr'll) = 1 (55) 

The integration goes over periodic commnting x^(r) and anti-commnting 'ip^ir), defined 
on the circle of the length T. 


Appendix B 

Let us consider the gauge invariant phase factor along small closed contour C for a 
particular inhomogeneous field 


W{C) 


<h(x, x) = Tr P exp 


A^dZfi 


The naive expansion over fields has the following form: 


( 66 ) 


z 

W{C) K. l + j) dz^ j> duyAi Ay(u)A^{z) + ... 
c c 


(57) 


The second term is gauge non- invariant (in nonabelian case). If one 
gauge ISni where 



j^dza{s) dz0{s) 
^ ds dx^ 




rewrites it in contour 


(58) 


the answer would be contour dependent. One can rearrange the series m to make it 
manifestly gauge-invariant: 


W{C) ^ l+i‘^ J da^p{z) j daua{u) Tr F^p(z, x^°^)Fj,^(m, x^°^) (59) 

s s. 


where = z{s = 0) and the shifted field strength is given by 

F^,(z,x(°)) = d>(xW,z)F^p(z)<I>(z,x(°)) (60) 

This is nothing but the nonabelian Stokes theorem |ddll^ . Notice that despite the Wilson 
loop itself is gauge-invariant and x’^°^-independent (it depends only on the contour C but 
not on the integration surface S in dSHD), the second term in the right hand side of (jsni) 
is gauge-invariant but contour dependent. Its SVZ-like expansion over condensate and 
derivatives would manifestly depend on the particular choice of contours Zp{s) and will 
be different ioT different contour gauges (i.e. for different choices of Za{s) and hence S). 
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However, in Gaussian approximation the account for contour dependence would not 
be legitimate ra. The reason is that the variation of the nonlocal Gaussian correlator 

over contours brings additional powers of Fafi-, and hence it is proportional to the correla¬ 
tors of higher orders (see related discussion in EP). We have used the same phenomenon 
replacing (EH) by (j21D in the main text. Indeed, let us look at the power expansions of 
those functions. They are given by 

!{s) = (TV (^1 + sDl + (61) 


and 


/.(i,) = {lVF, 




1 -h sD^ + — 


— -DaiFapDiJ + g-^a/3 


+ 


F,.) 


(62) 


It is seen that the difference between /(s) and / 2 (s) contains correlators of higher orders, 
resulting from non-commutativity of T)^’s. In Gaussian approximation each power of 
derivative corresponds to large factor and in this sense to extract the leading term 
one can always neglect the commutators [DqD/j] . The validity of Gaussian approximation, 
in its turn, is controlled by the dimensionless parameter 

K={TTF^^yX^ 

which is assumed to be small (it was mentioned above that, for example, iA(0)A^ is 
about 1/10 according to the lattice data). In physical terms, Gaussian vacuum is short- 
correlated, and the nonperturbative gluon condensate is small in units of the lightest 
gluelump mass (the latter is of the order of 1/A). Using the language of covariant pertur¬ 
bation theory one can say that we have summed up only terms of the kind {Tt FF) 
and not the terms with lesser number of derivatives. 


To summarize, we have 

/(s) = /2(s) + 0 (k‘+'') 

where r is some positive number. Therefore the leading Gaussian asymptotic of f{s) is 
given by the asymptotic of / 2 (s) found in eq. ()d2|l in the main text. 
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